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In this letter, atom optic techniques are proposed to con- 
trol the excitation of a Bose-Einstein condensate in an atomic 
trap. We show that by employing the dipole potential induced 
by four highly detuned travelling- wave laser beams with ap- 
propriate phases and frequencies, one can coherently excite a 
trapped Bose-Einstein condensate composed of ultracold al- 
kali atoms into a state rotating around the trap center. The 
connection to vortex states is discussed. 



The experimental realizations of Bose-Einstein conden- 
sition in an atomic trap |l) have led to a broad interest 
in the properties of coherent ultracold atomic gases. The 
comparison of trapped condensates to superfluid Helium 
suggests the possibility to study and analyse the rota- 
tional properties of the condensate This is usually 
made by exciting states in which the atoms move collec- 
tively around one or several vortex lines. Recently the 
properties of vortex states of a trapped Bose-Einstein 
condensate has been studied theoretically ||] . 

In this letter we propose to employ the atom optic 
techniques to study the excitation of rotating states, be- 
ing a coherent superposition of vortex states, from the 
ground state of the trap. In atom optics, the center- 
of-mass motion of atoms can be manipulated by light- 
induced gradient forces. We apply the same idea to 
a trapped Bose-Einstein condensate which can be de- 
scribed by a single macroscopic wave function. The inter- 
action of light waves with ultracold atomic samples com- 
posed of Bose-Einstein condensates has well been studied 
in recent years [Q- 0]. The vector nonlinear stochas- 
tic Schrodinger equations describe the dynamics of the 
ground- and excited-state atomic wave packets in the 
laser fields 1^. For large laser detuning, adiabatic elim- 
ination of the excited-state component leads to a non- 
linear scalar Schrodinger equation for the ground-state 
macroscopic wave function ipg of the condensate 
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The parameter A is connected to the s-wave scattering 
length of the atoms in the presence of the laser field. 
For simplicity we ignore the effects of the nonlinearity 
by setting A = in this paper. The incorporation of 
the nonlinear interaction in the proposed scheme will be 



the subject of a further publication. The center-of-mass 
Hamiltonian Hem. includes the kinetic energy and the 
harmonic potential of the trap. According to the usual 
treatment of the quantum mechanical harmonic oscillator 
the Hamiltonian can be expressed as 

Hem. = Tiuj^{al.ax + a^ay + 1} + hu)z{a\az + 1/2} , 

(2) 

where := [xi i = x,y, z, is the 

annihilation operators for the harmonic oscillator in the 
i-direction and Ri :— ^fij (Muji) is the trap size pa- 
rameter. M is the mass of the atoms, lu±_ is the trap 
frequency for the x- and y-direction, and lOz is the trap 
frequency in the z-direction. Because the trap is cylindri- 
cally symmetric around the z-axis, the Hamiltonian Hem. 
commutes with the orbital angular momentum operator 
J3 — xpy — ypx- This is a natural result of the conserva- 
tion of orbital angular momentum along the z-directon 
in the cylindrically symmetric trap in the absence of the 
laser fields. To excite the rotating states from a Bose- 
Einstein condensate occupying the ground state of the 
trap which is the ground state |0) of the harmonic oscil- 
lators, one must introduce an additional external force 
in the x-y plane to break down the rotational symme- 
try of the ground state around the z-axis and create an 
angular momentum for a rotation of the ground state. 
The laser- induced dipole potential Vl {x) in equation (|l|) 
can be used for this purpose. Before working out the ex- 
act form of the laser-induced dipole potential required for 
the excitation of vortex states, we need to understand the 
properties of the rotating states for the ground-state con- 
densate. A vortex state rotating around the axis of sym- 
metry is a state in which the macroscopic wave function 
ijjgix) has the form ex.p{inip)x{r, z), where ip is the angle 
and r is the radius in the x-y-plane and n is an integer. 
It therefore can be characterized as being an eigenstate 
of J3 with eigenvalue nh, and of Hem. as well. To deal 
with rotating states it is convenient to use the complete 
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set of states \v) = {cLY {c\_)™\0) /-s/ZIm! to decomp 
the wavefunction of the condensate. Here |0) denotes the 
common ground state of the three-dimentional harmonic 
trap, and the vortex annihilation operators c± are given 

by 
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(3) 



Using their commutators with Hem. and J3 one can show 
that c± both annihilate one quantum huj±_ of phonon en- 
ergy and change the angular momentum by ±?i. 
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In terms of the properties of the rotating states, we 
conclude that if the laser-induced dipole potential con- 
tains the appropriate combination of the vortex opera- 
tors (c±,c!|_), one can create a superposition of vortex 
states from the ground-state condensate. To obtain such 
a dipole potential, we propose to employ four travelling- 
wave laser beams with the spatial configuration as shown 
in Fig.l. 

The dipole potential induced by these laser beams has 
the following form 

4(A + z7/2)^^^ 

where we assume that all lasers are highly detuned from 
the transition frequency between the internal ground and 
excited state. In this case the excited state has only a 
small population and can be adiabatically eliminated. In 
equation (Q) the symbols 0.a{x) := d- Ei^\x)/h , a — 
1, . . . , 4 denote the Rabi frequencies for the laser beams 
with frequencies Ua- d is the matrix element of dipole 
moment of the two-level atoms, and Ea^'' (x) is the posi- 
tive frequency part of the ath laser's electric field, and the 
average detuning of the lasers is given by A. We assume 
uja — oJb to be much smaller than A so that the difference 
between the true detuning of a laser and A is of higher 
order in 1/A. The spontaneous emission rate 7 describes 
the incoherent loss of atoms from the condensate. The 
loss can be neglected in the time scale considered here 
for the case A 7. 

For the configuration shown in Fig.l, the two laser 
beams propagating along the x- and y-direction have the 
same frequency wi = ^2 = w. Their phases are chosen 
to be equal at the origin so that the Rabi frequencies of 
the two lasers have the form fli(x) — ^oexp[ikx] and 
^2(2^) = SIq exp[ifc?/]. The other two laser beams are 
chosen to have the frequency 0^3 = CJ4 = to' and prop- 
agate in directions slightly different from the x- and y- 
direction with a small angle 9. The frequency difference 
Alli :— uo — uj' is assumed to be so small that the dif- 
ference in the wavelength between the two pairs of laser 
beams can be neglected. In terms of the geometry shown 
in Fig.l, the wavevector of the third laser beam is given 
by kcos[d)ex — ksm(6)ey w ke^ — 6key, where Sk ;= kO 
denotes the deviation of the wavevector. Similiarly the 
4th laser's wavevector is given by key — Ske^- The phases 
are chosen so that laser 4 is in phase with lasers 1 and 
2 at the origin whereas the phase of laser 3 is shifted by 
— 7r/2. This phase shift is essential for the scheme to work 
since in the final Hamiltonian it will lead to a 7r/2 phase 
difference between and Oj,, which exactly realizes the 
vortex opertors (^. The third Rabi frequency is then 
given by ^■i{x) — exp[— i7r/2]r2o exp[i(fca:: — Sky)] and the 
fourth by ^4^{x) — JIq exp[i(fc?/— (5fca;)]. This arrangement 
of laser beams produces the dipole potential 
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j^e'^i'-^-y^+'^''''] } -f H.C. . (5) 

To deduce explicitly from Eq. an interaction that 
can create rotating states we consider the case that the 
laser intensity is so weak that the effective Rabi frequency 
|r2oP/(4A) is much smaller than the trap frequency ijJj_. 
This results in the following consequences: (i) If the 
atoms absorb and emit photons within a pair of laser 
beams with the same frequency (first square bracket in 
Eq. (||)) transitions between states of different energies 
Hem. are off-resonant and can therefore be neglected; (ii) 
If the atoms absorb a photon from one pair of laser beams 
and emit it into the other pair with different frequency, 
we have to take into account the frequency difference 
Aoj between the two photons with different frequencies. 
It produces an energy shift so that transitions between 
states having the energy difference ±?iAix' are resonantly 
enhanced whereas other off-resonant transitions are sup- 
pressed. Here we choose 

Aw = wj_ (6) 

so that only transitions between neighboring states of the 
harmonic oscillators are dominant. 

Mathematically this can be incorporated by switch- 
ing to the interaction picture with respect to -ffc.m., 
so that the annihilation operators are replaced by 
ai exp[— iwj^i] for i — x,y. In the interaction picture, 
by performing the rotating wave approximation with re- 
spect to Lu±, the exponential exp[ik{x — y)] has the form 

i 

e'''<'''-y^ « 6-"'/^ : Jo(2r/y^) : (7) 
e-*^'"*e''=("~^) « le^"'/' : Ji{2r]^/^) : , (8) 

where for notational convenience we have introduced the 
operator q := {ax^o-y)/ V^- The normally ordered Bessel 
functions are defined as 

V 2 / ^-^ mlim + n ! 

^ m— 

(9) 

Since the initial state for the problem discussed here is 
the condensate occupying the ground state of the trap, 
the Lamb-Dicke parameter 77 :— kR± is generally larger 
than one for the current realizable condensates |Q| which 
usually have transverse sizes R± > 1/zm. As a re- 
sult, Eqs. ^ and (||) imply that the terms containing 
eicp[±ik(x — y)] in Vl of Eq. are exponentially sup- 
pressed for large i] and therefore can be omitted for real- 
istic cases. A physical interpretation of this effect can be 
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given by realizing that exp{ikx) = exp{ir]{ax + aj,)) acts 
as a displacement operator on the trap states. Thus, for 
a large value of 77, it would cause a large displacement of 
the initial ground state of the trap. But since the inter- 
action is weak the energy required for this displacement 
is not fully provided so that only the low-energy part of 
the displaced ground state can be realized. 

On the other hand, the Lamb-Dicke parameter Sr] := 
SkR± corresponding to the exponentials exp[iSkx] and 
exp[iSky] in Eq. depends on the angle between the 
laser beams through Sk = kO and can therefore be ad- 
justed over a wide range. Here we require 9 to be very 
small so that (5?7 <C 1 is valid. This allows us to approxi- 
mate the two exponentials by 1 -I- i5kx and 1 -I- i5ky^ re- 
spectively. The small angle between the laser beams lead 
to a simple dipole potential Vl which linearly depends on 
the combination of the creation and annihilation vortex 
operators of the trap. 

In addition, in the interaction picture, 
the term exp[— z(5u;t]{l -I- i5kx} can be approximated by 
iSkRj_al./\^ under the rotating-wave approximation. A 
similar expression is obtained for exp[— zfJcj^Jjl + iSky} 
so that the final expression for the complete Hamiltonian 
in the interaction picture reads 
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(10) 



The physical content of Eq. (|T^) is very clear. The sec- 
ond term in Eq. ( |l0| ) exactly excites atoms in the ground 
state |0) into the rotating c_ mode. Therefore the initial 
cylindrical symmetry is broken by the interaction Hamil- 
tonian ( [l0| ) and a rotating state of the center-of-mass 
motion of the ground state atoms in the trap is created 
with each atom gaining an angular momentum h Thus, if 
we start from the ground state |0) of the trap, this term 
produces a vortex state with angular dependence exp(i(y9) 
in the coordinate representation. The term proportional 
to c_ is the corresponding Hermitian conjugate term and 
decreases the angular momentum by h. By ignoring the 
unimportant constant phase shift in Eq. (p^, and by 
returning to the Schrodinger picture, the time evolution 
operator U (t) of the condensate has the simple form 

{/(i) = exp |^[cL exp(zcj^t) — c_ exp(— icj^t)] j- , (11) 



where the time factor 
T„ := 
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characterises the time scale required for the creation of 
rotating states. When the time t is much smaller than 
Ty, the ground state evolves into the state U{t)\0) which 
approximately is a linear superposition of the first ex- 
cited vortex state and the initial condensate occupying 



the ground state |0) of the trap. With increasing interac- 
tion time t, higher vortex states are created. Finally we 
have a coherent state for the vortex operator c_ for a long 
interaction time. This rotating coherent state is a coher- 
ent superposition of vortex states circulating in the same 
direction at different angular velocities. From the coor- 
dinate representation of the vortex states, (a?|(cL)"|0) = 
2"/2 {x + iy)"i?;[^"+^^ exp[-(a;2 + y^)l{2R\)] , we can de- 
duce the probability density for the coherent rotating 
state, 
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p(t) = |V^(f|t/(0|0)P = ^cxp{ 



{x~XQ{t)f 



(13) 



where x^it) = (x) = {R±t/Ty)[cos{Lu±t)ex — sin(wj_i)ey 
determines the center-of-mass trajectory of the conden- 
sate. Eq. ( |l^ ) describes a condensate starting at the 
center of the trap and rotating with increasing amplitude 
thereby preserving its shape (see Fig. 2). 

So far we have theoretically presented a simple scheme 
to create rotating states from a trapped ground-state 
Bose-Einstein condensate by a light-induced dipole po- 
tential. To analyze the experimental feasibility, we use 
the experimental data of ^^Na Bose-Einstein condensate 
created by the MIT group ||l^ as an example. The mass 
of the ^^Na atoms is M — 3.8-10^^^ g and the wavelength 
of the ^■^Na D-line transition is 589 nm so that the wave 
number of the laser beams is given by A: = 1.06 • lO^cm. 
For the MIT trap the trap frequency in the x-y-plane is 
LUj_ w 27r • 320Hz so that the transversal trap size parame- 
ter is found to be R± ~ 1.17/im. With these parameters 
the Lamb-Dicke parameter 77 — kR± is approximately 
12.4 so that the right hand side of Eqs. (|^,H) can safely 
be ignored for low lying trap states. Hence the ^^Na 
Bose-Einstein condensate is an appropriate candidate to 
satisfy the conditions required above. 

The small Lamb-Dicke parameter Sr] = k6R± is deter- 
mined by the angle 9 and we can realize a value ^77 = 0.1 
by adjusting 9 around 0.46 degrees. In addition, to make 
sure that only resonant interaction is included, the effec- 
tive Rabi frequency ^eff '■= l^oP/(4A) must be much 
smaller than the trap frequency U!± . This can be experi- 
mentally realized if the laser intensity is chosen so weak 
that ^eff is about 200 Hz. In this case, the characteristic 
time scale for vortex excitation Ty is about 50 ms which is 
experimentally reasonable since it is much shorter than 
the life-time of the currently realized trapped conden- 
sates. 

In conclusion we have shown that it is possible to cre- 
ate rotating states from a trapped ground-state Bose- 
Einstein condensate in a trap by using light-induced 
forces. Technically we show that dipole potential induced 
by four travelling-wave laser beams with an appropriate 
configuration in space, phase and frequency can be used 
to realize such a purpose. The important features of the 
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proposed scheme are the following: (i) The low inten- 
sity of the laser beams to guarantee that only resonant 
transitions can occur; (ii) The small frequency difference 
between the two pairs of lasers to provides the energy dif- 
ference required to bring the atoms from the ground state 
of the trap into a superposition of the first excited states, 
(iii) The small angle 9 between the laser beams avoids 
the problem that the Lamb-Dicke parameter 77 is so large 
that the transitions are strongly suppressed if the differ- 
ence between the absorbed and re-emitted momentum is 
not very small, (iv) The phase difference between the 
laser beams leads to a special form of the light-induced 
potential Vl which exactly results in a creation of vortex 
states. 

Although the nonlinear interatomic interaction be- 
tween the atoms is ignored in our discussion, it can be 
shown that the total Hamiltonian including the nonlin- 
ear interatomic interaction commutes with the second 
quantized form of the angular momentum operator J^. 
Hence in principle the interatomic interaction does not 
affect the creation of rotating states, and only affects the 
spatial profile of the macroscopic wave function in the 
radical direction. The detailed discussions on the effect 
of the nonlinear interatomic interaction will be covered 
in a long publication ||^. 
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Figure Captions: 



Fig. 1: Proposed scheme to create a rotating state be- 
ing a coherent superposition of vortex states of the con- 
densate. The laser beams 1 and 2 are oriented along the 
X- and y-axis. A frequency modulator produces a slight 
frequency difference between the beams (1,2) (solid line) 
and (3,4) (dashed line). A small angle between beam 1 
and 3, and beam 2 and 4 is arranged by the adjustable 
mirror. 



Fig. 2: The trajectory of the center-of-mass of the con- 
densate in the x-y-planc under the influence of the four 
laser beams described above, x and y are given in units 
of the trap size R±. During the rotation the shape of 
the condensate remains a Gaussian profile with a width 
of 2RjL. 
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